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We calculate the collisional energy loss suffered by a heavy (charm) quark created at a finite time 
within a Quark Gluon Plasma (QGP) in the classical linear response formalism as in Peigne et al. 

We pay close attention to the problem of formulating a suitable current and the isolation of 
binding and radiative energy loss effects. We find that unrealistic large binding effects arising in 
previous formulations must be subtracted. The finite time correction is shown to be important only 
for very short length scales on the order of a Debye length. The overall energy loss is similar in 
magnitude to the energy loss suffered by a charge created in the asymptotic past. This result has 
significant implications for the relative contribution to energy loss from collisional and radiative 
sources and has important ramifications for the "single electron puzzle" at RHIC. 
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I. INTRODUCTION 

One of the experimental signatures of the creation of 
novel QCD matter at the Relativistic Heavy Ion Collider 
(RHIC) is the detailed suppression pattern of high trans- 
verse momentum hadrons This "jet quenching" has 
been previously explained as caused by the energy loss 
due to induced gluon bremmstrahlung in a hot Quark 
Gluon Plasma (QGP) 0. Recent data for the suppres- 
sion of single non photonic electrons [|[ (thought to probe 
the energy loss of heavy quarks), however, cannot be re- 
produced via purely radiative energy loss arguments with 
a physically reasonable set of parameters [J]. 

A recent publication by Mustafa @ stated that a rea- 
sonable set of parameters at RHIC would lead to elastic 
energy loss of a heavy quark that is of the same order as 
radiative energy loss. This conclusion goes against pre- 
vious estimates of the collisional energy loss which used 
asymptotic arguments to assert that the radiative energy 
loss is much larger than elastic @. Inspired by this de- 
velopment, a recent paper by Wicks et al. Q gets better 
agreement with the single non photonic electron data by 
incoherently adding the energy loss due to elastic as well 
as radiative sources while also averaging over geometric 
and gluon number fluctuations. It is extremely impor- 
tant to take fluctuations into account as they allow for 
the simultaneous reproduction of both the light quark 
suppression and the non photonic electron data. Their 
paper, however, neglects the finite time creation effects 
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on the collisional part of the energy loss that arise from 
the fact that the heavy quark probe is part of a back 
to back hard scattering occurring at a finite time in the 
medium. 

Peigne et al. [1| were the first to study the effects of 
the finite creation time of the probe on the collisional 
energy loss suffered by it. They claimed that the colli- 
sional energy loss was largely suppressed in comparison to 
the infinite limit for medium lengths applicable to RHIC 
physics. In the current paper, we redo the collisional en- 
ergy loss calculation in the same formulation as Peigne et 
al., paying close attention to the exclusion of all radiative 
effects and especially binding effects that arise if the full 
color current is approximated incorrectly. 

Recent papers by Djordjevic '9] and Wang [l(| ad- 
dress this problem using a quantum mechanical formal- 
ism rather than the linear response approach used in the 
current paper. The effect is modeled by restricting the 
interaction Hamiltonian to a finite time range in Djordje- 
vic 1 9] while Wang considers the interference between 
radiative and elastic amplitudes in first order medium 
scattering. The results from these two papers are in- 
consistent, however, since the full medium propagators 
suppress interference between radiative and elastic am- 
plitudes [9j. Our results are consistent with Djordjevic 

1. 

In Section |TT] we describe the subtle problems with for- 
mulating and interpreting currents that are created at a 
finite time. In Section HTT1 we formulate and derive calcu- 
lations for the infinite time case so as to show agreement 
with previous calculations and illustrate the subtleties in 
our specific method. In Section IIVI we reproduce the an- 
swer derived by Peigne et. al. and then correct it for 
radiative and binding effects. In Section [V] we present a 
summary of our results. 
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II. A COMMENT ON CONSERVED CURRENTS 

When one calculates elastic energy loss in the formal- 
ism of a classical macroscopic field theory, the general 
method is to posit a source current which travels in the 
medium in question and induces an electric field. The 
induced electric field acts back onto the current and does 
work on it, leading to a decrease in the energy of the 
source particle. 

In our application to jet production in nuclear col- 
lisions involves modeling a 2 — > 2 collision process in 
pQCD, with one of the final state particles being de- 
tected as a hadron after losing energy in the medium. A 
classical current that corresponds to this is, 

Jphys = ©(*) vxt) + q%v$6(x v 2 t b)) 

+0(-t) (g?<<5(x - vi't) + <7>/<5(x - va't - b)) . (1) 

The current in Eq. [T] has two particles with color charge 
ql and (Defined via q^q^ = flyC^a^) that experience 
a collision at t = with impact parameter b. Cr ; is the 
quadratic Casimir for representation Ri (Cr ; = |(3) for 
quarks (gluons)), a s is the QCD coupling constant (set 
to a s = 0.3) while wf = dox? — (1, Vi) is the four velocity 
of the i th particle. The impact parameter is determined 
by the momentum transfer in the hard collision Q, the 
'primed' velocities are the pre collision velocities of the 
particles, the 'unprimcd' velocities are the post collision 
velocities of the particles consistent with energy and mo- 
mentum conservation, and the final color configuration 
for the particles is the same as the initial configuration 
due to the Abelian approximation in force throughout 
this paper. 

We proceed with the analysis by calculating the in- 
duced field due to Eq. Q] and act with it on the particle 
whose energy loss is in question (say particle 1). The 
induced field will include creation radiation and interac- 
tions between particle 2 and particle 1 alongside the effect 
that we hope to calculate: the induced field of particle 
1 acting back onto particle 1. The radiative energy loss 
correction in QCD was first investigated in [llj and will 
be removed from consideration in this paper. The effect 
of particle 2 on particle 1, what we call the binding or 
shared interaction, averages out to zero in the following 
way. In order to calculate a single particle observable one 
has to sum over all the color orientations of particle 2 for 
each orientation of particle 1. This summation cancels 
the effects of the interaction and gives a zero effect on 
average. 

Due to the complicated nature of the current in Eq. [TJ 
early calculations [l3j in the induced field formalism ig- 
nored the effects of radiation, binding and finite time of 
the hard interactions and used a simple infinite time sin- 
gle particle current. 

^ (1) (x,t) = sW(x-v*) (2) 

The subscripts in Eq.[5]specify that this is an infinite time 
current that was created in the asymptotic past and that 



there is only one particle present. The current can also 
be written in Fourier space, 

j^ (1) ( W ,k) = 2^V^-k-v). (3) 

One can immediately see that the continuity equation 
written in Fourier space is satisfied, i.e. k^j^ ,^ = 
where k^ — (u>, k). 

In order to reintroduce finite time creation effects into 
the problem, we write down the naive generalization of 
Eq.H 

$J(x,t) = qV6(x-vt)e(t) 

$U«,k) = ( r V " - . ) . (4) 

K > \u> — k ■ v + in J 

Note, however, that the one particle finite time current 
is explicitly not conserved, as ^'^o k^j^ = iq a . This is 
expected because creation of charge out of the vacuum is 
a violation of charge conservation. 

The simplest generalization to a conserved finite time 
current is, 

i$(x,i) = g a «<5(x - Vl t) - <£(x - v a *))0(t). (5) 

This current models a current neutral pair production at 
t = 0. Two oppositely charged particles are created at 
t = and propagate outwards from the origin at pos- 
sibly distinct velocities Vi and V2. The Fourier space 
expression for this model current is, 

J(2)( w . k )=^ r — r — —J 

y 1 \u) — k • vi + n] uj — k • V2 + ir\ J 

The two particle finite time current satisfies charge con- 
servation and one can explicitly see that ^!^q k^j^. = 0. 
The values of the velocity cancel out in the calculation of 
the charge conservation condition so we can always set 
Vx = v and v 2 = 0. This changes the current to the 
following: 

J( ° 2 >,k) = iff J— 1-) 

iq a k • v 

lu + in uj — k • v + it] ' 



Eq. [7] was used as the current by Peigne et al. Q . 

Note that in formulating a suitable finite time current 
we were forced by charge conservation to include a second 
particle along with all the complications of the binding 
interaction. In order to study purely the finite time ef- 
fects, we will now have to remove the unrealistic binding 
terms as well as radiative effects. Recall that in the ac- 
tual physical situation corresponding to Eq.[T]the binding 
effect averages out to because the two particles have un- 
corrected charge in scattering of color neutral hadrons. 
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In Eq. [7] , however, the creation of the second particle 
out of the vacuum leads to exact anti correlation of its 
charge with the first particle, giving a large unrealistic 
binding. This effect is further enhanced by the fact that 
the second particle is created as a stationary charge at 
the origin which means that its fields are not relativis- 
tically contracted and its effects are longer range and 
stronger. We will show below that the enhanced unre- 
alistic binding to the second particle is the cause of the 
large suppression effect claimed in Peigne et al. Q. 

We will use Eq. as our current when doing the finite 
time analysis, taking great care to differentiate between 
self interactions of particles and shared two particle in- 
teractions. We also define an infinite time two particle 
current to further illustrate the effect of the interaction 
between multiple particles while avoiding the subtleties 
of the finite time problem. 



i£,( 2 )(",k) 

JL,(2)(^ k ) 



2irq a {8{uj - k • v) - S(lo)) 
27rg a v5(a;-k-v) (8) 



III. COLLISIONAL ENERGY LOSS WITH 
INFINITE TIME CURRENTS 

One can write down the rate of change of energy suf- 
fered by a particle traveling in a linear medium due to 
the induced field as, 



dE 
dx 



^e-^)^v.E»,(,k) (9) 



It is understood that v — |v|, // = k • v/kv and L is 
the path length that the particle has traversed in the 
medium. In order to regulate classical divergences (due 
to the Coulomb logarithm), the integration over Fourier 
modes will only be carried out to a maximum value k max . 
We assume k max is the maximum momentum transfer in 
a 2 — > 2 collision process between the incident particle 
and a medium parton in the center of mass frame, 



k m ax = min { E, 



2q(E+p) 



\JM 2 + 2q(E + p) 



(10) 



E is the initial energy of the incident particle, M is its 
mass, p is its momentum and q is the typical momentum 
of a medium particle, all in the medium (plasma) frame. 
For numerical calculations in this paper, we will consis- 
tently use M — 1.2 GeV in order to be representative of 
the charm quark energy loss and set q ~ T = 250 MeV. 

The induced electric field is the total medium field mi- 
nus the vacuum field produced by the current, E° nd = 
E a — E° ac . We can use the calculations in Appendix [5] 
to write the following expression for the produced fields. 



e q - ^( W ,k)(i-kk)j a + k£»|j( w ,k)j 00 

= D R {w,k)(wf -kf a ) (11) 



The propagators Dr{u, k), D^(w,k) and Z?^(w,k) are 
the vacuum Maxwell propagator, the in medium trans- 
verse propagator and the in medium longitudinal prop- 
agator respectively, all with retarded pole prescriptions 
and defined in Appendix [X] in terms of the Hard Ther- 
mal Loop (HTL) polarization functions, ILi^k) and 

n ± (w,k). 



A. Infinite Time One Particle Current 

We first calculate the energy loss in the infinite time 
one particle case in order to compare with previous calcu- 
lations [1, HI, EH ■ This means using Eq. [2] as the source 
current in Eq. 1111 Instead of getting an explicit expres- 
sion for E" nd , we separately calculate the total medium 
energy loss and vacuum energy loss felt by the particle 
and take the difference at the end. 

The current dependence of the vacuum field can be 
written as wjj* — ,^ — 2irq a (ujv — k)S(ui — k • v). 
Using the expression for Dr(uj, k) from Appendix lAl 



dE v 



dx 



.C F a s 
i (1 



v 2 ) 



dkdfj, ^ 
v z fi z - 



= 0. 



(12) 



The angular integration leads to the vacuum energy loss 
being equal to zero. 

We can express the electric field in Eq. [Tj] as the fol- 
lowing. 



E>, k) = \D^u, k)-kk + Djifa k)(l - kk) j 

(13) 

We use the expressions in Appendix [X] for the longitu- 
dinal and transverse medium retarded propagators and 
have used continuity of the current in order to write the 
th component of the current in terms of the longitudinal 
current. This expression for the field can be entered into 
the formula for energy loss. 



dE 
dx 



' v2n 2 



dk- 



k 2 



1 _ k 2 v 2 - (k- v) 2 
e II k 2 — (k • v) 2 ej_ 



(14) 

eu and ej^ are the longitudinal and transverse dielec- 
tric function respectively and are defined in Appendix 
lAl Both the dielectric functions are evaluated at uj = 
k • v + irj. This expression can be rewritten as the fol- 
lowing. 



dE 
dx 


CfO-s f „ k • v ., ... , . . 
v2 ^ k2 (W + /u(k)) 


h||(k) = 


/m(eii) 

hi 2 


h±(k) = 


{k 2 v 2 - (k- v) 2 )(k- v) 2 /m(e_L) 


\{k 2 - (k-v) 2 e ± | 2 



(15) 



The parts of the integral that depend explicitly on the 
real parts of the dielectric function vanish by symmetry. 
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FIG. 1: (Color Online) The figure shows the energy loss frac- 



tion 



ae _ r 



dE/dx 



as a function of the incident momen- 



tum p. Curves are shown for an incident charm quark with 
M = 1.2 GeV and L — 5 fm. It includes the full calculation as 
well as previous approximations from the literature; Bjorken 
[f|, Thoma- Gyulassy [13 ] and Braaten-Thoma [13 ] . 



In the leading logarithm approximation (as in [13|) we 
neglect the effects of the Debye screening mass in the 
denominator of the integrand in Eq. [15] and evaluate the 
integral analytically. 



dE 
d.r 

m 



C F a, 



, log 



m D 



1 (v -1(1 -v 2 ) log 



(16) 



Here 



47ra s T 2 (l + N F /6) is the squared Debye 



screening mass. We use Np — for calculational pur- 
poses which is equivalent to assuming that the plasma is 
purely gluonic. Notice that the infrared singularity that 
was regulated by the full dielectric functions in Eq. [15] 
is apparent in Eq. [16] and has been regulated "by hand" 
at the scale mp. The ultraviolet regulator k rnax is the 
same as before. The logarithmic dependence on energy 
(via kmax) is characteristic of collisional energy loss and 
is maintained even in the full calculation when one does 
not make the leading logarithm approximation. 

The full calculation can be compared to the leading log 
approximations to the collisional energy loss that have 
been used previously in the literature. 



dE A 
dx 



CfC( s 



m 2 D f(v)log(B A ) 



(17) 



The index A runs over the different calculations at- 
tempted. They all agree with each other up to the kine- 
matic ratio inside the logarithm Ba = ^ ma3: , given in Q • 
We denote them as Bjorken (Bj) Q, Thoma-Gyulassy 



(TG) Q and Braaten-Thoma (BT) V, 
B Bj = (AET)/{m 2 D ) 
4pT 



£?tg — 



B 



BT 



(E- P + 4T)> /imD) 
2^7 0.85 ET^j I (2^) 
4^7 0.36 ^)/ (4) E«Ml 



(18) 

For purposes of calculation, the turnover in Bbt is set at 
E = 2.36-^-; this value is determined by the condition 
that Bbt should be a continuous function of E. 

Figure Q] shows the energy loss fraction — as a func- 
tion of the initial momentum of the incident particle with 
the medium length traversed L — 5 fm. Curves are shown 
for the full calculation of dE/dx as in Eq. [TSland for all 
the previous leading log formulae for collisional energy 
loss detailed in Eqs. [T7] and [15] We include curves for 
the full calculation at two different values of k ma x ■ There 
is some ambiguity in setting the value of the momentum 
of the medium particle q ~ T — ST. Both curves are 
shown as a measure of the maximal theoretical uncer- 
tainty. Note that all the calculations agree for asymp- 
totically large energies (as they should) and that the full 
calculation (Eq. [T5"|) does not possess the unphysical en- 
ergy gain at low momentum that is shown by all the 
leading log approximations. 

The agreement between the full calculation in this pa- 
per and the earlier leading log calculations is a non triv- 
ial result. In particular, the Bjorken and Braaten-Thoma 
calculations are attempted in an entirely different formal- 
ism than our own. Their formalism calculates collisional 
energy loss as an average over the quantum mechanical 
cross section for elastic collisions in the medium. The 
formalism of the current paper and in Thoma-Gyulassy, 
however, uses linear response arguments in the abelian- 
ized classical Yang-Mills theory. 

The agreement between our calculation and Thoma- 
Gyulassy is also a subtle one. One can explicitly see 
by comparing Eqs. Ql)] and [T7] that our calculation is a 
factor of 2 larger than Thoma-Gyulassy, modulo different 
kinematic ratios in the logarithm. This factor of two (also 
noticed by Peigne et al. [3) comes from the different 
polarization functions n_j_(o;,k) and Hm (c*j, k) (defined in 
Appendix [XJ used in this paper as compared to the ones 
used in Thoma-Gyulassy. Thoma-Gyulassy use only the 
positive frequency half (u> > 0) of the logarithmic branch 
cut present in the polarization function while the current 
paper takes into account the full branch cut extending 
across —k<ut<k. Thus, the Thoma-Gyulassy result is 
suppressed by a factor of 2. 

Their result is enhanced by the use of a large k max in 
the kinematic logarithm as compared to all other calcula- 
tions. They evaluate the maximum momentum transfer 
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in an elastic collision in the plasma frame rather than in 
the center of mass frame, thus including a 7 factor as 
an enhancement in their result. This enhancement can 
easily be seen by calculating the asymptotic energy de- 
pendence of k m ax in each model. The Thoma-Gyulassy 
kinematic factor for asymptotic energies gives k max ~ E 
while all others give k max ~ y/ET. We also note the cal- 
culation by Djordjevic Q uses the same kinematics and 
branch cut characteristics as Thoma-Gyulassy and thus 
the momentum dependence of energy loss in the Djordje- 
vic model is very similar to the Thoma-Gyulassy model, 
despite the fact that in spirit and calculational formalism 
Djordjevic is closer to Bjorkcn or Braaten-Thoma. 



B. Infinite Time Two Particle Current 

We now employ the two particle infinite current as in 
Eq. [5] in order to clarify the effects of binding without the 
complication of a finite time analysis. The end result of 
the finite time analysis will remove all reference to this 
binding effect. We do the calculations in this section on 
our way to make contact with the calculation of Peigne 
et al. 8], showing how including this unrealistically en- 
hanced contribution leads to a suppression in the energy 
loss. 

The current in Eq.[5]includes components of the energy 
loss that can be thought of as a self interaction (equiv- 
alent to the one particle case) while the rest is a shared 
interaction between the stationary second particle and 
the moving particle in question. We will call terms of 
the first type self terms while the second kind will be 
shared terms. 

In order to facilitate easy physical interpretation of 
the different sources of the energy loss, it is conve- 
nient to decompose the retarded propagator in terms of 
a purely advanced propagator Da, and a difference of 
the retarded and advanced propagator Z3_, i.e. Dr = 
Dr — Da + Da = D- + Da- The procedure is detailed 
in Appendix O This helps to isolate the singularities in 
the complex lu plane that are present in the propagator 
and thus manifestly breaks up the expression in Eq. [5] 
into contributions from different sources. This is particu- 
larly useful for the finite time current to be studied later 
and will be used in this section for purely illustrative 
purposes. 

We will calculate the vacuum and medium field energy 
losses separately, paying close attention to which parts 
of the energy loss are coming from true feedback in the 
plasma and which are due to the interactions with the 
stationary particle. 



1. Vacuum Field Contributions 

We start by looking at terms in the energy loss that 
come from the vacuum field specified in Eq. 1 1 1 L The 
expression for the retarded propagator Dr{uj, k) can be 



taken from Appendix IA1 and its decomposition into Dr = 
D- +Da is detailed in Appendix [Cl The energy loss from 
fields due to the _D_(w,k) propagator are as follows. 



dE- 



dudk 



dx J (2tt) 4 
£>_(u;,k)(a,j^ (2) (a;,k) 



]y£ ( 2)(«,k)) 



(19) 



As explained earlier, this contribution can be broken up 
into a self interaction and a shared interaction. Using 
Eq. H we write wj^ >(2) - kj£ (2) = 27rg°((wv - k)<y(w - 
k ■ v) + kS(ui)). The part proportional to 5(u> — k ■ v) 
contributes to the self interaction while everything else is 
a shared interaction. We separate them, 



dE-^ f C R a 



dx 



J rp— ,shar 

dx 



^(i-, 2 )|fc^( M -£) 



5 ( [J, H — ) }dkdfi 







= -C R a s J k 2 fie tk ^ L (d(-k) - S(k))dkdn 
= (20) 



As shown in Appendix [Cl the vacuum D_ contribution 
is just twice the radiative pole contribution. For an in- 
finite time current, the only kind of radiation possible is 
Cerenkov radiation. The characteristic Cerenkov cone is 
present in the self interaction expression as /1 = cos(#) = 
~. This contribution goes to zero as the velocity of the 
particle is less than that of light. 

The only contribution left over is that due to the ad- 
vanced propagator DA(w,k). This can be written as, 

dE vac _ <f f d^dk -i<LU-kvu)± „ 



dx v J (2tt) 4 
£> A (u;,k)( W v • j" (2) (w,k) - v ■ kf" (2) (u,k)) (21 



The self interaction part integrates to zero but there is a 
non zero shared interaction. 



dE. 



Adv,shar 



dx 



CfOL s 



d 3 k 



v .ke ikllL D A (Q,k) (22) 



This vacuum shared interaction is a measure of the 
binding energy of the pair of particles and is the only 
source of vacuum energy loss for the moving particle. 
It can be thought of as the energy lost by the moving 
particle as it climbs out of the binding potential well of 
the stationary particle of opposite charge. This integral 
is nominally infinite as both particles appear at the ori- 
gin at t = 0. The divergence is simply due to the in- 
finite shared energy of two point charges that goes like 



kmax ~ ■ It is typical of the binding (shared) terms 

seen in the analysis and should be regulated using the 
length scale |b| (and hence momentum scale Q) of the 
initial hard interaction. In order to be consistent with 
our simplified current which has no information about 
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an initial hard interaction and with the medium contri- 
butions that have their own prescribed interaction scale, 
we will use the previously defined medium k max as our 
regulator. 



2. Medium Field Contributions 

The medium field can be conveniently broken up into a 
transverse and a longitudinal contribution and those in- 
dividually can be further subdivided into a D_(cj, k) con- 
tribution and a Da(lu, k). As shown in Appendix [Cl this 
decomposition leads to energy loss contributions from 
three different physical processes: the particle-hole cut, 
the plasmon radiation loss, and the advanced coulombic 
field. 

We begin by considering the longitudinal electric field 
in the medium (as given by Eq. As in the vacuum 
case, the only kind of plasmon radiation possible is of 



Cerenkov type, which is not present as > k (as 

shown in Appendix [C|l . We therefore move on to contri- 
butions from the particle- hole branch cut. 



dE cut,self 

dx 



C F a s 



d 3 k 

(2^)3 



v-L4H(k-v,k) (23) 



The function A" (u), k) is 
part of Dl(cj,k) that is 
particle- hole branch cut. 
frequency for — k < u> < 
the particle hole cut intej; 

There also remains an 
causes energy loss. This 



defined in Appendix [Cl It is the 
linked back to the logarithmic 
Thus, it only has support in 
k. The shared contribution to 
;rates to zero since A" (0, k) = 0. 
advanced longitudinal field that 
contribution can be written as, 



dE 



Adv,self 

i 



dE 



dx 

Adv .shar 



= c 



dx 



-CFOt & 



(2tt 
d 3 k 



d 3 \< - ii 
^v.kDi(k-v,k) 



(27T) 



3 v-ke tkvi ^(0,k).(24) 



d\(u),\l) is the longitudinal in medium propagator with 
advanced boundary conditions. 

Due to the form of the current ( 2 )( w : k), the trans- 
verse field only provides self interactions and no shared 
interactions. As in the longitudinal medium and the full 
vacuum cases, the radiative part of the field has to be 
purely Cerenkov in nature and thus does not contribute. 
We begin with the transverse particle- hole branch cut, 



dE cut,self 

dx 



= C 



F a s v 



d 3 k 

(2^)3 



(l-(v-k) 2 )A ± (k-v,k) (25) 
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FIG. 2: (Color Online) The figure shows the differential en- 
ergy loss —dE/dx as a function of the path length traversed 
in the medium L. Curves are shown for an incident charm 
quark with M = 1.2 GeV and initial energy E — 10 GeV. 




FIG. 3: (Color Online) The figure shows the energy loss —AE 
as a function of the path length traversed in the medium L. 
Curves are shown for an incident charm quark with M = 1.2 
GeV and initial energy E = 10 GeV. 



The only remaining piece is the advanced transverse field. 

-C F a s W^4(l-(v-k) 2 )£>i(k-v,k) (26) 



dx 



(2tt) 3 



D^(w,k) is the transverse in medium propagator with 
advanced boundary conditions. 



A ± (u, k) is defined in Appendix [Cl and is the part of the 
D^(ui, k) propagator that can be associated with the log- 
arithmic branch cut singularity structure of the particle- 
hole continuum. Similar to the longitudinal case, this 
function only has support in frequency for —k<uj<k. 



3. Total dE/dx 

We now add up all the energy loss due to the medium 
fields and subtract all the vacuum energy loss expressions 
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FIG. 4: (Color Online) The figure shows the differential en- 
ergy loss —dE/dx due to the binding between the two particles 
that constitute the incident current as a function of the path 
length traversed in the medium L. Curves are shown for an 
incident charm quark with M = 1.2 GeV and initial energy 
E = 10 GeV. 



In order to make sure that this depression in the en- 
ergy loss is coming purely from the difference of shared 
interaction we can systematically remove the effect. 



dE, 



(2) 



dx 



E dELd 



dx 



dE a 



(1) 



dx 



E 

AG{self} 



dEt 



dx 



(28) 



If we leave out only the vacuum interactions but keep 
all the two particle interactions in the medium we get 
dE(2)/dx shown as the red dashed curve in Figs. [2] and 
[3] This result makes sense as it is just the one parti- 
cle infinite result with the additional energy loss due to 
binding between the two particles in the medium. We 
attain dE^/dx by removing all mention of the shared 
interaction entirely. This result goes back to the infinite 
time one particle it should. This exercise is a 

good proof in principle that we can perform the neces- 
sary subtractions in the finite time case in order to isolate 
the effect we want to study. 



in order to determine the total work done by the induced 
field. 

dEmd _ dE^ led _ dE„ ac 
dx ^— ' dx ^— ' dx 

A B 

The indices A, B run over the different types of energy 
loss (e.g. from the branch cut, from advanced fields), 
both self and shared. Figs. [2] and [3] show this energy 
loss. The black dot-dashed curve in the figures shows 
the differential energy loss (dEi n d/dx) as calculated in 
Eq. [27l The blue solid curve is the full calculation for 
the infinite time one particle energy loss as shown earlier 
in the paper. Note that the induced field energy loss (in 
Fig. \2§ is actually smaller for some time than the one 
particle energy loss. 

The initial dip in dEi n d/dx (seen in both Figs. [2] and 
[3]) is caused by the subtraction of the vacuum contribu- 
tion from the medium parts of the energy loss as stated 
in the induced field prescription. Why this should lead 
to a dip in energy loss is seen in Fig. [U Here we plot 
only the shared (binding) terms that contribute to the 
problem. Due to screening effects, the in medium energy 
loss is less than the vacuum energy loss. In our calcula- 
tion for dEi n d I dx in Eq. [27] we have taken a difference 
between these binding energies causing a small energy 
gain; a "Ter-Mikaelian" adjustment to the binding en- 
ergy. This is exactly the effect that is included in Peigne 
et al. ; an energy gain from the difference of binding 
energies between vacuum and medium interactions. The 
magnitude of this component is also enhanced because of 
the unrealistic charge correlation between the two parti- 
cles in the current and by the fact that the second particle 
is created at rest at the origin. 



IV. THE FINITE TIME PROBLEM 

The subtleties involved in choosing a conserved cur- 
rent for the finite time case are manifest in this section, 
which means that we must start with the two particle 
finite time current j^ 2 "(w,k). We will analyze the full 
energy loss, specifying exactly how to remove the effects 
of the second particle which cause the anomalously small 
results of Peigne et. al. Q as well as the radiative parts 
that are included in our energy loss calculation (as shown 
in Appendix |B|) . Note that we will not combine the two 
separate terms inside j^(w, k) into one term but rather 
will treat each term distinctly as a self and shared inter- 
action. 



A. Induced Field Energy Loss and Removing 
Binding Effects 

1. Vacuum Field Contributions 

We begin the analysis by considering the contribution 
from the advanced propagator D^Wjk). The analysis 
is greatly simplified by the fact that all the singularities 
in the propagator are in the upper half of the complex 
lu plane while the current singularities all reside in the 
lower half. We can thus perform a contour integral in 
frequency space and close the contour in the lower half 
plane, picking up the residues at the current singulari- 
ties. This calculation leads to the same advanced field 
energy loss as in the infinite time case (Eq. l2Tj) . as the 
current pole has the same frequency dependence as the 
infinite time current. This will be true of all the advanced 
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contributions in the paper. 



dE- 



Adv. shar 



dx 



c 



FCt s 



d 3 k 

(2tt) 3 



Mce' fcflL £» A (0,k) (29) 



We can now move on to the D_{uj,k) contribution. 
The self and shared contributions are, 



dE~' self 



.■ °f as I Am+ 



d J k(r+ +r s 



r 
r 



sei/ fc jfc - k • v 

e i(fe+to M )i fci;2 + v . k 



self 



k k + k • v 



(30) 



=-^/^(r+ ar + r^ r ) 



^(2tt) 2 

r 



e -i(fc-fcu/j)i v , k 



r 



e i(k+kvfi) v v • k 



shar 



k k 



(31) 



Unlike the infinite case, there will actually be a non zero 
contribution from the D_(w,k) propagator. This is eas- 
ily interpreted as the radiation emitted when a charge is 
created out of the vacuum. 



2. Medium Field Contributions 

The contribution from the longitudinal particle hole 
cut is, 



dE cut,self k A ll (a;) k)v.k 

— — ^ — — ZO^C^s / — — — 



dx 

cut, shar 



(27T) 



— k • v + i?7 



gator. It can be written down as 



dx 



pii.-t- 



— i(uA-fcuu) — i 



1 self ~ 



dE pl,shar 



c^-k-v 
v • k 



/(-<>k) o^+k-v 



(34) 



P i 



dx 



4-+ 



shar 



v(2tt) 2 J ^ shar shar) 

—Uu",—kvu,)— i 

e v p' ' " v • k 

/(^ ; ,k) ^jr 

e «(Wp(+^)# v . k 



/(-^/,k) w 



(35) 



pi 



These plasmon contribution are those field components 
that can support the longitudinal plasmon dispersion 
relation relation wjj ; (fc), plotted in Appendix [B] The 
function /(w,k) is the longitudinal plasmon wave func- 
tion renormalization factor and is defined in Appendix 
|B1 in terms of the derivative of the polarization function 

n„(w,k). 

There are no shared contributions from the transverse 
fields, only self interactions. Contributions from the 
transverse particle-hole cut can be written as, 

dE cut,self „ duidk _ i(u _u vlA h 



dx 



(2tt) 



x A^,k)(l-(v.k)») (36) 
lu — k • v + irj 

We can also get the advanced contribution very easily 
from contour integration. 



dE 



Adv,self 



dx 



fdudk i(u _ k ^ AH(q;,k)v ■ k 

J (27T) 4 UJ + 17] 



dx 



Cfol s v 



d 3 k 

J2nf 



(l-(v.k)^)£>i(k.v,k)(37) 



The integrals in the previous expressions can be peformed 
by doing the angular integrals first in terms of the Expo- 
nential Integral function. The D^(co,k) contribution is 
easy to calculate via contour integration in the uj plane. 



dE* dv - self r d 3 k 

C F a s I ^v-kflUk.^k) 



The plasmon terms coming from the transverse parts of 
the field can be written as, 



dE 



pi, self 



dx 



dE 



dx 

Adv,shar 



dx 



= -a 



(2^ 

d 3 k 



-uL,+ 
self 



self 



^ 2 -^± ( rf3 k( r^+ + r^-) 

v{2tt) 2 J V self selt ' 
w^e-*^-* ")- l-(vk) 2 



k- v 



^e'^+fc^)- l_( v . k ) 2 



(38) 



There is now a non zero contribution in the finite time 
case from the plasmon term inside the D_(ui, k) propa- 



The function g(u>, k) is the transverse plasmon wave func- 
tion normalization factor and is defined in Appendix [B] 
in terms of the derivative of the polarization function 

n ± (w,k). 
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FIG. 5: (Color Online) The figure shows the differential en- 
ergy loss —dE/dx as a function of the path length traversed 
in the medium L. Curves are shown for an incident charm 
quark with M — 1.2 GeV and initial energy E — 10 GeV and 
include radiative modes. 
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FIG. 6: (Color Online) The figure shows the energy loss —AE 
as a function of the path length traversed in the medium L. 
Curves are shown for an incident charm quark with M = 1.2 
GeV and initial energy £7 = 10 GeV and include radiative 
modes. 



in Figs. [5] and [6] The solid blue curve in the figures is 
the answer to the infinite time single particle problem, 
presented for comparison. The total induced field energy 
loss (counting both self and shared interactions) is shown 
as the dotted line. Note that the induced field energy loss 
is far less than the infinite case and does not get back up 
to the infinite answer (in Fig. [6]) even after traveling u 
to 10 fm. This calculation reproduces Peigne et al. 
and is different only in as much as the parameter sets 
used are different. 

As explained in the previous section, a naive calcula- 
tion of energy loss from the total induced electric field 
includes shared interactions as well as the self ones. The 
shared contribution is unrealistic and needs to be system- 
atically removed from the calculation. We first remove 
the vacuum shared field, 



dE, 



(2) 



dx 



A 



A 

med 



dx 



E 

Be{self} 



dEi 



dx 



(40) 



We refer to the energy loss that neglects the binding dif- 
ference as dE(2)/dx as it explicitly shows the effect of 
the stationary particle's in medium potential well on the 
moving particle. The dashed curves in the Figs. [5] and 
IH1 show this situation, which is just the one particle fi- 
nite contribution plus the energy loss due to in medium 
binding. 

In order to compare the finite time calculation with 
the equivalent infinite time answer, we need to remove 
all the shared parts (medium as well as vacuum). 



dE a 



(i) 



dx 



E 

Ae{self} 



dx 



E 

Bd{self} 



dEi 



dx 



(41) 



The result of the purely self interactions is shown as the 
dot-dash curves in Figs. [5] and [51 The finite time cal- 
culation starts lower than the infinite time case as the 
induced field takes some time to build up and be created 
in the medium. We can now see in Fig. [5] that the dot 
dash-curve (which still includes radiative energy loss ef- 
fects both gluonic and plasmonic) is getting closer to the 
infinite time case. We have systematically removed the 
unphysical binding terms and shown that they were the 
cause of the large suppression of energy loss in Peigne et 



3. Total 

The first object we calculate is the energy loss due 
to the total induced field, making no distinction as to 
whether the field is coming from a self or shared interac- 
tion. 



dE, 



ind 



dEi 



dEi 



dx 



A 



dx 



B dX 



(39) 



The indices A and B sum over all the different sources of 
energy loss, both self and shared. The results are shown 



B. Subtracting the Radiative Ter-Mikaelian Effect 

We note that the one particle result quoted in the pre- 
vious section includes more effects than just the pure 
collisional energy loss suffered by the particle. When a 
current is created out of the vacuum, it carries with itself 
certain off shell Fourier modes that can be radiated off 
to infinity in the form of gluons in the vacuum and plas- 
mons in the medium. As shown in Appendix[Bl the Joule 
heating expression we use includes contributions that be- 
come the total zeroth order in opacity radiative energy 
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FIG. 7: (Color Online) The figure shows the differential en- 
ergy loss —dE/dx as a function of the distance traveled L 
due to radiative effects at zeroth order in opacity. Curves are 
shown for an incident charm quark with M = 1.2 GeV and 
initial energy E = 10 GeV. 
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FIG. 8: (Color Online) The figure shows the energy loss —AE 
as a function of the distance traveled L due to radiative effects 
at zeroth order in opacity. Curves are shown for an incident 
charm quark with M = 1.2 GeV and different initial energies 
E = 10,15,20 GeV. 



loss when we integrate dE/dx for an infinite path length. 
This radiative energy loss needs to be removed from the 
calculation if we are to correctly compare the finite time 
result with previous estimates for the collisional energy 
loss. 

The radiative effects manifest themselves in terms that 



dE;. 



dx 



id] i 



can be schematically written down as 
where the index A can represent a vacuum, medium 
transverse, or medium longitudinal contribution and the 
function D rad (u>, k) is the radiative Green function as 
specified in Appendix|SJ The radiative energy loss curves 



for a particle with initial energy E — 10 GeV in both vac- 
uum (solid blue curve) and medium (dashed and dash- 
dotted curves) are shown in Fig. [7] There is only trans- 
verse gluon radiation in the vacuum but the medium sup- 
ports both transverse and longitudinal polarizations for 
the plasmons. The longitudinal energy loss, however, is 
significantly smaller than the energy loss caused by trans- 
verse plasmons. Note that the vacuum peak is larger than 
the transverse peak. 

The Ter-Mikaelian effect (first studied by Djordjevic 
et al. [11]) is the difference between the zeroth order 
radiative energy loss in the medium and the zeroth or- 
der vacuum radiative energy loss (shown in Fig. [7]). We 
implicitly include this difference in the calculations for 
dE(x)/dx at the end of the last section when we calculate 
the induced field (E a — E° ac ). It is shown in Appen- 
dices [B] and [C] that the radiative part can be isolated by 
carefully dealing with the singularities in the complex uj 
plane. The result for the radiative part is, 



dE 



Bad 



dx 



K A£{pl} 



dE ted 

dx 



dE, 



B 

vac 



se{-} 



dx 



(42) 



The factor of a half comes from the Z?_ (u>, k) propagator. 
This difference propagator has a part that has the correct 
radiative dispersion relation but has a residue of twice the 
radiative residue in the retarded propagator. 

The difference between medium and vacuum radiative 
energy loss at E = 10 GeV is shown as the solid black 
curve in Fig. [7J Note that this curve includes the ef- 
fects of the longitudinal radiation (something which is ne- 
glected in Djordjevic et al. [ll|) and also has a complete 
calculation of the wave function renormalization factor 
(Appendix [BJ . Fig. [S] shows the total Ter-Mikaelian en- 
ergy radiated by the particle as a function of the path 
length L for different values of the initial energy. One 
can see that this effect is small (~ 2% ). It does, how- 
ever, have a positive correlation with the initial energy 
of the particle. 



C. The One Particle Analysis 



The finite time energy loss can be schematically rep- 
resented as AE ~ {DR e t — Dji a d)j. The removal of the 
radiative poles explicitly subtracts all radiative contribu- 
tions. We further write the retarded propagator in terms 
of the advanced and mixed parts, AE ~ (DAdv + D- — 
Dnad)j- As shown in Appendix [Cl the mixed propagator 
carries all the pole structure of the retarded propaga- 
tor, only with twice the residue. So we can further write 
AE ~ (D Adv + A + 2D Rad - D Rad )j. The function A is 
the cut contribution specified in Appendix [Cj Using this 
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FIG. 9: (Color Online) The figures show the differential energy loss —dE/dx suffered by a particle traveling in a medium as 
a function of the path length traversed L for different values of the initial energy E. The solid blue line is the infinite time 
calculation, the red dashed and dash-dotted curves are the stable and transient field energy losses respectively and the black 
dashed curve is the total finite time energy calculation. Radiative modes are excluded from consideration. 




FIG. 10: (Color Online) The figures show the energy loss —AE suffered by a particle traveling in a medium as a function of 
the path length traversed L for different values of the initial energy E. The solid blue line is the infinite time calculation, the 
red dashed and dash-dotted curves are the stable and transient field energy losses respectively and the black dashed curve is 
the total finite time energy calculation. Radiative modes are excluded from consideration. 



procedure we can write the total finite time answer. 



dx 
dx 

dEgfable 

dx 



dx 



dx 



2 \ -V dx 

K Ae{ P i} 

, dx 

A^{cut,adv} 



E dEyg C 
r , dX 

Be{-} 



(43) 



We have broken up the one particle differential energy 
loss into a 'transient' part and a 'stable' part. The 'tran- 
sient' field consists of the single set of simple poles that 
have been left over after the removal of the radiative 



energy loss. The 'stable' fields consist of the contribu- 
tions coming from the particle-hole cut and the advanced 
part of the propagators. Both these curves are shown in 
Figs. 151 and [TU1 along with the total energy loss and the 
infinite time calculation for comparison. The curves are 
shown as a function of the path length L and for different 
values of the initial energy E. 

The 'transient' fields are called as such because they 
are created at small initial times of the propagation in 
the medium and then quickly die away (as seen in Figs. [3] 
and [TO]). They tend to have a large effect since the cur- 
rent in use is created instantaneously out of the vacuum 
(with 'infinite' acceleration). Djordjevic [9( finds similar 
contributions in the quantum mechanical calculation and 
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attributes them to the lack of time translation invariance, 
another way of invoking initial conditions. 

The 'stable' field starts off at zero and then steadily 
builds up to the infinite case, reaching about 90% of the 
infinite field strength by a time on the order of one De- 
bye length as one would expect (Fig. O . It is this delay 
in the set up of the fields by a time on the order of the 
Debye length which is the real retardation effect. Fig. ITD1 
shows the integrated energy loss as a function of the path 
length traversed for different initial energies. As one can 
see, the total energy loss in the finite case is of the same 
order as the infinite case and even slightly enhanced over 
the infinite case. This enhancement comes from the con- 
tribution of the 'transient' field effects that have not had 
time to settle down. The 'stable' field energy loss has a 
small depiction initially but builds up after a path length 
on the order of one Debye length and then the slope of 
the curve becomes equal to the infinite case. This result 
is in agreement with the complementary calculation per- 
formed by Djordjevic 0. A calculation with the current 
as formulated in Eq. [T]would be useful in order to achieve 
an estimate for the magnitude of the transient effect. 



CONCLUSIONS 



We find that close attention must be paid to the form 
of the current used in the analysis using linear response 
formalism. Some finite time currents lead to problems 
in interpretation and implementation of the formalism 
due to subtleties involved in ensuring that the current is 
conserved and that all radiative and binding effects have 
been discarded from the comparison being made to infi- 
nite time results. The unphysically large binding effects 
are particularly important for the form of the current 
used by Peigne et. al f§[ , as a failure to account for them 
can cause a spurious interpretation of the results as a 
large time delay in the onset of collisional energy loss. 

Our results show that there exists a small time delay 
(on the order of the Debye length) during which the chro- 
moelectric fields build up to their full intensity. Thus 
there is an initial depletion in energy loss in the finite 
time case as compared to the infinite case. 'Transient' 
field effects, however, push the energy loss further to- 
wards the infinite case; an effect that is sensitive to the 
detailed form of the current. These 'transient' fields, 
however, are unrealistically enhanced due to the form 
of the current and should be studied for more realistic 
charge configurations in order to estimate their true ef- 
fects. Our results are consistent with the calculations 
of Djordjevic 0], providing a complementary approach 
that both verifies and provides further physical insight 
into the result for non asymptotic heavy quark energy 
loss. 
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APPENDIX A: SOLUTIONS TO MAXWELL'S 
EQUATIONS 

The solution to classical Yang-Mills theory in the 
Abelian approximation is equivalent to solving the 
Maxwell equations of electromagnetism. In the vacuum 
the inhomogeneous part of the Maxwell equation can be 
written as the following. 



47rf 



(Al) 



Here, F^ a = d»A va - d u A^ a is the field strength tensor 
and A^ a is the four vector potential (gluon field). Note 
that the non linear part of the field strength has been 
omitted due to the Abelian approximation. Equation 
IA1I can be solved in the Lorentz gauge and in fourier 
space to yield A^ a (cj,k) — ~iDji(uj,k)j^ a (uj,k) where 
we use the retarded Maxwell propagator, Dn(u),k) = 
7 aT 4 y2 1 1,2 ■ The electric field can be further calculated 
in Fourier space as the following. 



E; oc ( w ,k)=^( w ,k)( w j»-k/ a ) 



(A2) 



In the medium one has to solve the Maxwell equa- 
tions with macroscopic fields in the linear response for- 
malism. We define the macroscopic fields D a = e -E a 
and B a =A* H a , with the permittivity and permeabil- 
ity tensors defined as the following: ^ (w,k) = 1 and 

e* (w,k) = e x (w,k)(l - kk) + e\\ (w, k)kk. We solve the 
macroscopic Maxwell equations in Fourier space for the 
electric field. 

E a (w,k) = Ei(w,k) + Ejf(w,k)k 

Ei(w,k) = D£(w,k)j5.(w,k) 

£(f(w,k) = £>f(a,,k)j >,k) (A3) 

The transverse and longitudinal propagators are defined 
in terms of the transverse and longitudinal permittivities. 



I>i(w,k) 
D L R {w,k) 



k 2 — (u> + irj) 2 e±(u> + if], k) 
—Airi 



fee||(w + irj, k) 



(A4) 



The permittivities in turn can themselves be written in 
terms of the HTL one loop calculation for the polariza- 
tion tensor, £_l(cg>, k) = 1 — IL±(lj, k)/u> 2 and e||(w,k) = 
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1 + II|| (u>, k)/k 2 . The transverse and longitudinal polar- 
izations arc. 



1 



n B (*,k) = ml, l--log — 



A; 2 /w 
log 



(A5) 



Here m 2 D = 4ira s T 2 (l + Np/6) is the squared Debye 
(screening) mass where T is the temperature of the 
medium and Np is the number of quark flavors that are 
in thermal equilibrium. 



APPENDIX B: JOULE HEATING AND 
RADIATIVE ENERGY LOSS 

Joule heating is defined as the work done on a current 
by the electric field at its local position. The total work 
done by the electric field on the current can be written 
as the following. 



W = J ctofxj Q (x,<) -E a (x,f) 



/ 



dudk 



r>,k)-E>,k) 



(Bl) 



This expression can be used to calculate the work done 
by an electric field on its own source. 

Specifically for the vacuum case we can use Eq. IA2I to 
write the following. 



W 



dudk 
(2^ 



D R (u,k)(u(\j ± \* + \jtf)-kjlf i ) ( B2 ) 



Note that all expressions quadratic in the current involve 
an implicit sum over the color index. We separate the 
retarded propagator into a pole residue and a principal 
value part in order to isolate the radiative parts of the 
field. 



D B (w,k) = PV 



D rad (uJ,k) 



4tt 2 



D rad (iv,k) = - — (6(u-k)-8(w + k)) (B3) 

We can isolate the radiative pole in the propagator in 
order to calculate the work done on the current due to 
the radiative field. 



dk 

dk 

4^2 1 



lin 



(B4) 



Eq. IB4I shows that the energy lost to Joule heating by 
the radiative part of the fields is exactly equal to the 
radiative energy loss of the current. Note that we need a 
conserved current in order to get rid of the longitudinal 
current parts in Eq. IB4I 
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FIG. 11: (Color Online) The figure shows the dispersion rela- 
tions for the radiative fields in the vacuum and medium. The 
calculation is done for T = 250 MeV, Nf = and a s = 0.3. 
Note that (0) 
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We can do a similar calculation for the medium case 
in which we have separate longitudinal and transverse 
contributions to the radiation. Using Eqs. IB 1 1 and IA3I 
we can derive the following. 



(27T) 



4^(^ k )ll 



wl{ = lw? Dku) > mf (B5) 



The medium retarded propagators can be separated into 
a radiative part and a principal value part in a way anal- 
ogous to Eq. 



Di(«,k) 

D. 



rod 



DUu,k) 



D" — - 



PV(Di(u,k))+Djr ad (u;,k) 
PV(D{(u;,k))+Dl ad (w,k) 

47T 2 fc 



(*(w-<4) + *(w + fcA))B6) 



Here the functions w^(fc) and uiK(k) are just the trans- 
verse and longitudinal plasmon dispersion relations re- 
spectively (shown in Fig. [TTj) . Note that the value of 
the plasmon frequency is always greater than the wave 
number (w^' 1 " > k). The functions g(u>,k) and f(u>,k) 
are defined as derivatives with respect to uj of the inverse 
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propagators. 
/(w,k) 



9 w (^_ fe 2_ n±( p 



2u 



2k 2 



1 / aH-fc 



9 w (fc 2 + n,|(-)) 



2 - /c 2 2fc 



1 , { u) + k 
log 



(B7) 



One can now determine how much work is done in joule 
heating due to fields of radiative type in the medium. We 
can show using Eq. IB6I inside Eq. IB5l that this amount is 
exactly equal to the zeroth order energy loss to plasmon 
radiation in the medium. 



W \, 

r r rad 

wL 



dk 2oj 
4^ 2 5 (c^,k) 
dk 2k 2 



|j±(^,k)| 2 



|j||(^,k)| 2 (B8) 



We have used the continuity equation in fourier space to 
represent the longitudinal radiative energy loss purely in 
terms of the longitudinal current. The expressions in Eq. 
IB8I include the wave function renormalization factors due 
to interactions in the medium as well as the dependence 
on the current. 



APPENDIX C: RETARDED, ADVANCED AND 
MIXED PROPAGATORS 

In order to maintain causal effects, all fields in the pa- 
per are calculated from the sources using retarded bound- 
ary conditions. Thus a field can be written schematically 
as E = Diij. We use a particular decomposition in order 
to separate the singularity structure of the propagator 
into something more amenable to physical interpretation. 
Lets first illustrate this in the vacuum case. We can write 
the retarded propagator in the following manner. 



D R (u,k) 
D_(w,k) 



D A (u,k) 

D A (uj,k) 



D A (oj,k) 



(CI) 



In the previous expression we introduce the mixed propa- 
gator (which is just twice the half retarded half advanced 
Feynman propagator) £)_(<j,k). The mixed propagator 
can be written as the following. 



D_(w,k) = D R (u,k) - D A (u,k) 
An 2 



k 



(S(uj -k)- S(uj + k)) (C2) 



Note that the advanced and retarded vacuum propaga- 
tors are only different when it comes to the position of 
singularities in the complex plane. Taking their differ- 
ence, as the mixed propagator does, just picks out this 
singularity structure. Note that this contribution is just 
twice the radiative contribution. We will see that this 
same procedure with the medium propagators will lead 
to the isolation of the particle hole cut as well as the 
plasmon pole. 

The basic difference between the vacuum and medium 
cases is the existence of the branch cut singularity in the 
medium propagators that is not present in the vacuum 
case. The principal value part that simply cancels out 
in the vacuum mixed propagator actually produces this 
electron-hole pair cut in the medium. 



A ± (u,k) 
A" (uj, k) 



A x (cj,k) + 2D^ d (u;,k) 

2Tr 2 m 2 D uj 2 (uj 2 - k 2 )e{k + u)Q(k - u) 
k 3 (k 2 -oj 2 + IL ± (^))(k 2 -uj 2 + n ± (^i)j 



A\u,k)+2Dl ad {u,k) 

~An 2 m 2 D ujQ{k + w)Q(k - uj) 

(fc 2 + n||(4H ))(fc 2 + n ll (^H)) 



(C3) 



The functions A _L (a;,k) and AH(w,k) occur as a result 
of the branch cut singularities present in the medium 
propagators. They can be interpreted as a contribution 
from the particle-hole continuum present in the plasma. 
Note the difference between these branch cuts and the 
cuts in the propagators used in Thoma-Gyulassy [l3| as 
only half the cut is used in that paper. 
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